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Метрология интеллектуальных систем 


В статье рассмотрены системные принципы метрологического обеспечения проектирования 
систем с элементами искусственного интеллекта. Структурно-стохастический принцип аппроксимации 
информации пространственно-временных полей обеспечивает преодоление существенной априорной 
неопределенности и высокие метрологические показатели по скорости адаптации, достоверности и 
устойчивости результатов. 


Введение 


При решении проблем, связанных с оценкой состояний широкого класса 
исследуемых объектов, с целью последующего принятия решений и выбора стратегии 
управления сложными процессами, все большее распространение получают методы, 
основанные на использовании информации полей различной физической природы. 
Значительное место среди этих полей занимают пространственно-временные поля, 
образуемые практически всеми объектами природной и техногенной среды в широком 
спектральном диапазоне. Путем использования информации пространственно-временных 
полей решаются такие задачи, как оценка качества теплоизоляции объектов с 
энергетической установкой, экологический мониторинг окружающей среды, медицинская 
диагностика и др. По мере расширения круга решаемых задач и их усложнения, к 
системам обработки информации полей предъявляются все более жесткие требования 
по основным метрологическим показателям их эффективности: высокая достоверность 
при ограниченном времени получения конечных результатов, устранение фактора 
субъективности при принятии решений, быстрая адаптируемость к новым задачам, 
высокая устойчивость результатов. 

Физические поля большого числа технических и медицинских объектов довольно 
разнообразны и сложны по своей структуре, а поэтому, как правило, требуют для 
интерпретации и последующего принятия решений специалистов высокой квалификации. 
Высококвалифицированные специалисты при анализе информации полей быстро 
утомляются, что в условиях дефицита времени приводит к недопустимому возрастанию 
числа ошибочных решений. Особенно проблематична безошибочная интерпретация 
оператором изменяющихся во времени реализаций полей, т.е. при проведении 
динамических наблюдений. Создание же автоматизированных систем, требующих 
предварительной настройки, связано с длительными по времени процессами изучения 
структурных особенностей полей и установления корреляционных связей между 
структурой полей и состояниями исследуемых объектов, т.е. с накоплением базы 
знаний. Как показывает практика, такой «ручной» способ накопления знаний крайне 
неэффективен и поэтому часто происходит отказ от решения той или иной 
практической задачи. 

В связи с этим актуальна проблема разработки принципов метрологического 
обеспечения адаптивного управления обработкой информации широкого класса 
пространственно-временных полей, обеспечивающего выполнение указанных требований 
в условиях существенной априорной неопределенности как о характере структуры 
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текущих реализаций полей и числе образуемых ими диагностических классов, так и 
об истинных состояниях исследуемых объектов. Суть данной проблемы сводится к 
созданию метрологического обеспечения методологии проектирования систем 
принципиального нового класса — адаптивных информационно-распознающих систем 
(АИРС), обеспечивающих возможность быстрого приобретения новых знаний о мало- 
изученных процессах. При этом поля, отражающие различные состояния исследу- 
емых объектов, предъявляются в режиме адаптации без указания принадлежности к 
тому или иному классу, а используемые в АИРС принципы управления обработкой 
должны обеспечить определение числа и характеристик образов полей в простран- 
стве их представления (описания) практически в реальном масштабе времени. 
Указанная проблема является составной частью общей, фундаментальной проблемы, 
связанной с созданием математических моделей систем с элементами искусственного 
интеллекта. Особую остроту данная проблема приобретает при обработке информации 
пространственно-временных полей, обладающих сильно выраженной изотропностью 
и образующих значительные области взаимного перекрытия диагностических 
классов (образов). 


Цель работы 


Разработка системных принципов метрологического обеспечения адаптивного 
управления обработкой информации пространственно-временных полей. 


Постановка задачи 


Идентификационные признаки информации реализаций полей можно 
представить в виде изменяющихся во времени совокупностей пространственно 
ориентированных градиентов, образующих связные конфигурации разных размеров 
(площадей). Такое определение признаков делает возможным реализацию принципа 
измеримости информации полей. Однако при этом существует целый ряд факторов, 
оказывающих мешающее воздействие: сочетание, неоднозначность и инверсный 
характер различных признаков, маскирующий эффект структуры полей. Действие 
перечисленных факторов и приводит к образованию значительных областей взаимного 
перекрытия диагностических классов (образов), что, в свою очередь, определяет 
вероятностный (стохастический) характер виртуального пространства функцио- 
нирования систем, а также принципов формализации процедур управления адаптацией 
и распознаванием образов полей. 

Если в качестве общей модели измерений использовать дифференциальное 
уравнение вида: 


В=Ф(В,5,0,1, (1) 


где В - вектор значений измеренных параметров идентификационных признаков 
(вектор сообщений); 5 - вектор состояний исследуемого объекта; И — погрешность 
измерений, то оптимальная оценка вектора 5, дающая решение задачи, будет 
определяться соотношением вида: 


^ 
$. = м. |, (2) 
В, — совокупность значений векторов сообщений в интервале времени [1:1]. 


Практическое применение соотношения (2) связано с поиском апостериорного 
распределения 5. Это очень сложная задача, решение которой в общем случае пока 
не найдено. Теоретически, согласно методам теории стохастических дифференциальных 
уравнений, для нахождения оптимальных оценок вектора состояний необходимо 
решить уравнение апостериорной плотности вероятности вектора состояния 5 после 
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получения результатов измерений. Решение о том, к какому из классов (образов) 
принадлежит данная реализация поля, обычно принимается по максимуму 
апостериорной вероятности. Однако следует отметить, что методов точного решения 
подобных уравнений пока не существует. Численное же их решение в практических 
задачах также оказывается невозможным из-за несуразно больших временных и 
вычислительных затрат. Применение приближенных методов ограничивается высо- 
ким порядком дифференциальных уравнений, при этом число уравнений быстро растёт 
с увеличением числа измеряемых параметров. Кроме того, в случаях стохастического 
(вероятностного) разделения полей на дискретное множество образов приходится стал- 
киваться с необходимостью формирования обучающих выборок несуразно больших 
размеров даже при относительно небольшом числе образов. Так, например, при общем 
(небольшом) числе алгоритмических констант, равном 20, необходимый размер 
обучающей выборки при разделении только на два образа (класса) окажется равным 
примерно 10°, чего уже вполне достаточно для отказа от использования стохасти- 
ческого подхода, а следовательно, и от решения указанной проблемы, тем более, что 
ситуация усугубляется отсутствием априорной информации как о реальном числе обра- 
зов, так и степени информативности признаков. 

Выход видится в использовании элементов структурного подхода, согласно 
которому всё множество пространственных градиентов делится по свойствам измери- 
мости и связности на непересекающиеся подмножества — элементарные составляющие, 
образующие ортогональный базис [1], [2] и являющиеся способом отражения (формали- 
зации) идентификационных признаков реализаций полей. Множество элементарных 
составляющих определяет, в свою очередь, допустимое множество образов полей. 
Каждая элементарная составляющая измерима через число её элементарных компонентов 
и число связей между ними. С помощью структурного подхода можно добиться суще- 
ственного снижения требований к объёмам обучающих выборок, т.к. каждой конкретной 
совокупности соответствует существенно меньшее число алгоритмических констант, чем 
их возможному множеству. Таким образом, решение поставленной научной проблемы 
возможно в рамках структурно-стохастического принципа аппроксимации виртуального 
пространства адаптивного управления обработкой информации полей. 


Метрология структурно-стохастического принципа 
управления 


Согласно данному принципу пространственно-временное поле можно рассматри- 
вать как некоторую случайную функцию следующего вида 


Ко = УДО) бьуи),  1=Ъ4, (3) 


где {0 (0)} — множество значений весовых коэффициентов 1-й элементарной составля- 
ющей; Я 7,1) — бинарное изображение 1-й элементарной составляющей; 4 — число 


составляющих. 

Данная модель позволяет реализовать многоуровневый иерархический метод 
структурного анализа, согласно которому происходит поочерёдное выделение составля- 
ющих с одновременным оцениванием их весовых коэффициентов и определением 
значений идентификационных параметров. В основу построения алгоритмов выделения 
элементарных составляющих положен принцип логической фильтрации. Результатом 
работы алгоритмов является принятие решения по каждому элементу поля о его 
принадлежности к той или иной составляющей. В основу построения алгоритма форми- 
рования изображений элементарных составляющих положен принцип объединения 
бинарных элементов в связные совокупности, который реализуется путем последователь- 
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ного наращивания множества без учета возможных направлений его распространения, 
т.е. применением алгоритма перебора в ширину. В основу определения значений иденти- 
фикационных параметров положен принцип измеримости информации полей. В итоге 
формируется вектор значений параметров (вектор сообщений). 

Обработка пространственно-временного поля производится на каждом такте {, ‚ т.е. 


в каждый момент времени поступления очередной реализации поля. По окончании 
обработки выдается результат в виде некоторого сообщения а, ={и;а,,....а,},, где 


п — число отдельных (парциональных) сообщений. Для всей последовательности 
реализаций поля может быть построена траектория изменения значений векторов 
сообщений во времени а,(#), представляющая, по существу, образ поля. В качестве 


образа структурно-подобного подмножества полей будет выступать не отдельная 
траектория, а некоторая их совокупность — трубка близких, в некотором смысле, тра- 
екторий А. Рассматривая каждую из полученных реализаций как некоторый случайный 


процесс, естественно использовать в качестве меры близости траекторий вероятностную 
меру — максимум апостериорной вероятности, обеспечивающий минимум средних потерь. 
Априорные вероятностные свойства любой траектории а, (+) могут быть заданы 


совместными конечномерными плотностями вероятности в точках ее отсчета 
в (и)ьараь) (а,,...а,)= те п (а... ах). (4) 


Апостериорную плотность вероятности траектории а (е ) представим как 


условную плотность распределения, когда уже получен ряд текущих векторов 
сообщений Ё,,...,Б, в моменты времени (,..., 1, ‚ в виде формулы Байеса 


7, =Р =Р(а,... а, В....В,) = 
Е Руль @зныа) 
А Бьзв.) 


где В(Ь,....Ь, |. а,) имеет смысл условной плотности вероятности распределения 


Е, (@1›--> Чу), (5) 


значений векторов В,...,Б,; В, (а,,...а,)— совместной плотности вероятности значе- 
ний векторов а,,....а,; Р(Б,....6,) — совместной плотности вероятности значений 
векторов Б,,...,Б, на /[-траектории. 

Совместная плотность вероятности значений векторов а,,....а, и В,...Б, может 


быть найдена путем интегрирования произведения в числителе соотношения (5) по 
а,....а; :. Знаменатель соотношения (5) представляет собой безусловную плотность 


вероятности значений фены ‚› КЮторая может быть найдена путем интегрирования того 
же произведения по всем а,,...,а,. Представляя условную плотность, входящую в 


числитель (5), на основании теоремы умножения вероятностей и выполняя указанные 
операции интегрирования, получим 


ВО = | В Боба вена ХВ, ланноа аа 
..хР(ЬЬ;а,а,)Р(Ь|а,)В(а,,...альа, )аа, аа 
| Оо аб нов В Бы Вань вы) ая 
..хР(ЬЬ;а,а,)Р(Ь, а, )Р(а,...а»а, аа, аа 


(6) 


«Штучний 1нтелект» 32008 21 


Анцыферов С.С. 
= 


Чтобы упростить последнее соотношение, будем полагать, что между сообщениями, 
получаемыми в последовательные моменты времени, существует вероятностная взаимо- 
связь, аппроксимируемая марковским процессом. В результате такой аппроксимации 
соотношение (6) принимает следующий вид 


Р(® = |. | ВФ а) вв, а, 1х... х Ва. )В а )Р(а)х 
хР(а, а, 1)х...х Р(а а, )аа...а,_ 
[Ро вор ид х-.х Рае) (а,)х 
хР(а, а, 1)х...х Р(а а, )аа....а, 


(7) 


После выполнения ряда преобразований, направленных на устранение операции 
многомерного интегрирования, получим 


Е, 8 
Е(ЮЕ (Ва, аа, й 


А 


& Р(ЮР(Ь 
п= ее 


где Р(®) имеет смысл априорной плотности вероятности принадлежности текущего 
сообщения траектории а, (+), а Ра) — функции правдоподобия. Результатом выпол- 
ненных аппроксимаций является стохастическое описание виртуального простран- 
ства — стохастические образы структуры векторов сообщений. 

Процедура распознавания является составной частью процесса адаптации. 
Случай перекрывающихся трубок (образов) только подчеркивает стохастический 
характер задачи, а предложенная стохастическая модель образа структуры векторов 
сообщений оказывается полностью ему адекватной. 

Для отыскания наилучшего значения вектора параметров образов р, миними- 


зирующего величину оценки среднего риска \(р,) и сходящегося при К -—> со к зна- 


* 
чению р, минимизирующего, в свою очередь, величину самого среднего риска, 
можно использовать рекуррентное соотношение следующего вида: 


Ре = Ре -Б, (6, ,а, У „(В рь в). (9) 


Здесь критерием оптимизации параметров образов служит градиент меры принадлеж- 
ности У „1, (В,, р, 1) к образу (8(В,,а,) — функция потерь). Данный критерий оптими- 


зации позволяет существенно ускорить процесс адаптации, особенно в случае 
перекрывающихся трубок траекторий, не меняя по существу режим распознавания, в 
котором главным критерием остаётся МАХ?) = МАХР . Важной новой составляющей 
процесса адаптации является формирование «следов» текущей траектории — результатов 
коррекции параметров существующих образов. Это позволяет путем определения 
меры принадлежности опорной траектории к «следу» и последующего определения 
градиента меры её принадлежности к опорному и скорректированному образам 
принять одно из решений: текущая траектория принадлежит трубке (образу) А, а ее 


след А! является скорректированным опорным образом 2 ; «след» текущей траектории 


является прообразом и в дальнейшем служит основой для формирования новой 
трубки траекторий (нового образа). 

Следствием взаимного перекрытия трубок траекторий (образов) является 
изменение уровня (меры) информативности векторов сообщений. Кроме того, в 
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процессе динамических наблюдений меры информативности элементарных составля- 
ющих поля, а соответственно и описывающих эти составляющие векторов сообщений, 
непрерывно изменяются за счет процесса усиления одних и ослабления других 
составляющих. В качестве оптимального правила оценки значений мер информатив- 
ности векторов сообщений, адекватного многокомпонентному характеру реализаций 
поля и учитывающего существование некоторого закона информативной связи 


между значениями парциального сообщения р, = |8" в различные моменты времени, 
можно использовать определение апостериорной плотности вероятности математи- 


ческого ожидания переменной идентификации ел в К-й момент времени при условии, 


что известны аналогичные оценки и значения сообщения В р во всех предыдущих тактах 


== Реаньв+.1) о 


Вводя марковскую аппроксимацию на последовательное изменение меры информатив- 
ности сообщения от одного момента времени до другого, получим взаимосвязанную 
пару рекуррентных соотношений, позволяющую находить оценки мер информативности 
как в данный момент времени, так и в следующий: 


56) Р(е! Ре! | р.) 


2 | 11 
[Ре Рем (1) 


где Р(е")= [ Ре" Р(е! |2 де. (12) 
Е 


Соотношения (11) и (12) обеспечивают возможность построения стохастических 
образов мер информативности каждой из элементарных составляющих и соответ- 


ствующих им векторов сообщений. Определение функции правдоподобия Р(2 *) 


осуществлялось исходя из того, что её вид должен соответствовать наихудшему виду 
функции распределения для каждого из сообщений (диагностических параметров), т.е. 
виду, приводящему к образованию зон взаимного перекрытия между трубками, к 
размытию границ между ними (рис. 1). Аппроксимация потока векторов сообщений 
марковским процессом определяет вид стохастического образа динамики этого 
потока. Наличие этого образа полностью раскрывает механизм селекции сообщений, 
что обеспечивает устойчивость режимов обработки информации. 

Оценка метрологических параметров структурно-стохастической модели АИРС 
проводилась как путем моделирования по методу Монте-Карло, так и путем обработки 
реализации полей реальных физических объектов. Как показали результаты 
моделирования, способ управления дискриминаторной характеристикой (градиентом 
меры принадлежности) с ограничением диапазона изменения межтактовой разности 
функции правдоподобия обеспечивает повышение скорости сходимости почти в 1,5 раза 
в случае двух соприкасающихся трубок траекторий и числе идентификационных 
параметров по каждой элементарной составляющей (число составляющих 5), равном 5. 
Радикальное улучшение условий адаптации в случае соприкасающихся трубок 
траекторий достигается при использовании способа управления с ограничением в 
сочетании с управляемым выбором корректирующих коэффициентов в рекуррентном 
соотношении (9). Исследование скорости сходимости при различных условиях взаимного 
перекрытия трубок траекторий показало, что ухудшение сходимости до некоторых 
пределов почти пропорционально увеличению степени перекрытия, после чего сходи- 
мость просто не достигается (кривая 4 на рис. 2). 
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Рисунок 1 — Функции распределения и правдоподобия в зоне перекрытия 
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Рисунок 2 — Зависимость сходимости от степени перекрытия трубок траекторий 
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Использование стохастического образа мер информативности для управления 
адаптацией в случае перекрывающихся трубок траекторий сводится к информативному 
«взвешиванию» стохастического образа структуры векторов сообщений путем 


умножения его на РЕ > | В таком случае 
Ут, к = Р(ЕЁ Р(ЮАР(а), (13) 


Тогда, с учетом влияния Р(®) и имеющего нормальное распределение Р(Е! ) с 


о, <в, <в;, дискриминаторные характеристики имеют вид, представленный на рис. 3. 


Введение пороговых ограничений по АР(Ба) предотвращает образование ложных 


следов траекторий, делающих процесс адаптации практически невозможным. 
Как показали результаты моделирования, учет меры информативности приводит к 
заметному ускорению сходимости в условиях взаимного перекрытия трубок траекторий, 
причем способ порогового ограничения векторов сообщений по мере информативности 
эффективен в широком диапазоне возможных уровней перекрытия. Проведенный 
анализ показал также, что для интенсивных потоков, содержащих большое число 
элементарных компонент, эффективно управление, связанное с селекцией векторов 
сообщений по информативности. Наличие в составе трубок неверно отнесенных к 
ним траекторий (погрешности адаптации) изменяет крутизну дискриминаторной 
характеристики, а следовательно, снижает скорость сходимости. 


— 57 
и |7 УПк. к-1 


Рисунок 3 — Дискриминаторные характеристики 


Вероятность ошибки распознавания с учетом погрешностей адаптации определяется 
следующим соотношением: 
1 


Б, =5(-рР--+РР-1 


а (14) 


где Р., - вероятность ошибки при отсутствии погрешности адаптации; 


ош 
Р - вероятность безошибочной адаптации. 
Данная зависимость представлена на рис. 4. 
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Рисунок 4 — Зависимость достоверности распознавания от погрешности адаптации 


Вводя в рассмотрение стохастический образ параметра семейства законов распреде- 
ления векторов сообщений и используя рекуррентное соотношение (9), можно найти 
оптимальные значения параметров данного образа для каждой трубки и для каждого 
такта траектории. В процессе моделирования данного режима было установлено улуч- 
шение сходимости в среднем примерно на 20 % и повышение устойчивости достовер- 
ности результатов (уменьшение вероятности ошибки) почти на 15 %. 

К наиболее реальным способам практической реализации АИРС можно отнести 
чисто программный и аппаратурно-программный с использованием микропроцессорных 
комплектов [3]. Эксперименты по использованию АЙРС для решения задач ранней 
диагностики опухолевых заболеваний молочной железы при проведении массовых 
профилактических обследований показали, что в результате адаптации обучающая 
выборка пространственно-временных тепловых полей оказалась разделенной на три 
существенные с точки зрения конечного результата трубки траекторий (рис. 5) мер 
принадлежности: трубка 1 — соответствует классу нормы, 2 — мастопатии, 3 — подозрению 
на рак. Использование результатов адаптации позволило почти на 20 % поднять достовер- 
ность термографической диагностики и приблизить ее к 90 %. 


Рисунок 5 — Трубки траекторий мер принадлежности 
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Вывод 


Предложенные системные принципы метрологического обеспечения структурно- 
стохастического управления обработкой информации пространственно-временных полей 
позволяют создавать адаптивные информационно-распознающие системы, способные 
функционировать в условиях значительной априорной неопределенности и пересекае- 
мости образов полей, обеспечивая тем самым высокие уровни достоверности и устой- 
чивости результатов, а также существенное ускорение процессов приобретения новых 
знаний о связях между состояниями исследуемых объектов и их физическими полями. 
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С.С. Анциферов 

Метролопя 1нтелектуальних систем 

У стати розглянут! системн! принципи метролойчного забезпечення проектування систем з елементами 
штучного 1нтелекту. Структурно-стохастичний принцип апроксимацй 1нформацй просторово-часових полв 
забезпечуе подолання 1стотно! апр1орно! невизначеност! 1 висок! метролопчн! показники швидкост! адаптаций, 
достоврност 1 стйкост! результат. 
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